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Abstract 

We measured the angle- and polarization-resolved reflection and photoluminescence spectra of 
ZnO inverse opals. Significant enhancement of spontaneous emission is observed. The enhanced 
emission not only has good directionality but also can be linearly polarized. A detailed theoret- 
ical analysis and numerical simulation reveal that such enhancement is caused by the abnormal 
dispersion of a high-order photonic band. The frozen mode at a stationary inflection point of a 
dispersion curve can strongly modify the intensity, directionality and polarization of spontaneous 
emission. 

PACS numbers: 42.70.Qs, 42.25.Fx, 78.55.-m, 78.40.-q 
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I. INTRODUCTION 



Photonic crystals (PhCs) have been extensively studied for control of spontaneous emis- 
sion. The first proposal by Yablonovitch in 1987^ utilized a complete photonic band 
gap (PBG) in a three-dimensional (3D) PhC, which was realized in 2004 by Noda and 
coworkers.- Compared to one- dimensional and two-dimensional PhCs, the 3D PhC has the 
utmost control of light-matter interactions. However, it is much more difficult to fabri- 
cate 3D PhCs, especially with the top-down approach. An alternative approach based on 
self-assembly attracts much attention due to ease of fabrication. Such a bottom-up ap- 
proach has been widely used for fabrication of face-centered-cubic (FCC) structures such 
as opals or inverse opals. An opal that is made of dielectric spheres has only partial PBG, 
namely, the PBG exists only within certain angle range instead of all angles. An inverse 
opal can have a full PBG if the refractive index contrast is high enough. It has been 
shown^'^'^i^'^'^'^ii^'ii'i^'^'^'^'^'^'i^'^>22i2ii22 that a partial PBG can effect the spontaneous 
emission process in an opal or inverse opal by suppressing spontaneous emission into certain 
directions. Enhancement of spontaneous emission is also observed at the edge of PBG. Most 
of these studies are focused on the lowest-order PBG in an opal or inverse opal. There are 
only a few experimental studies of the effect of high-order photonic bands on spontaneous 
emission.-i^i^ Recent studies2^i^'^i^'^i^i^'2^i2^i2ii2^i^iILi^i2£dS show that a high-order pho- 
tonic band can have abnormal dispersion which leads to many interesting phenomena such 
as super prism^i^i^ and negative refractive index— i^. In this paper, we demonstrate that 
a high-order photonic band with abnormal dispersion can significantly modify the photolu- 
minescence (PL) intensity, directionality and polarization in an inverse opal. 

The dispersion of photons in a PhC can be dramatically different from that in free space. 
Let us denote the dispersion of a photonic mode by uj{k), where u is the photon frequency 
and k is the wave vector. A mode with dtu/d/c ~ is called a slow mode because the group 
velocity is nearly zero. There have been many proposals of utilizing the slow modes of PhC 
to reduce the speed of light by orders of magnitude. A serious problem that hinders the 
slow light application is that a typical slow mode with (Puo/dk"^ ^ has large impedance 
mismatch at the PhC/air interface, thus the conversion efficiency of incident light into the 
slow mode is very low. To solve this problem, Figotin and Vitebskiy proposed to use the 
photonic mode at the stationary inflection point of a dispersion curve of a photonic band.— 
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Such a mode has both duj/dk ~ and (fuj/dk"^ ~ 0. It is called a frozen mode. When the 
incident light is in resonance with a frozen mode, the vanishing group velocity is offset by 
the diverging electromagnetic energy density. The energy flux inside the PhC is finite and 
comparable to the incident flux. Hence, the incident light can be completely converted to 
the frozen mode instead of being reflected. Our aim is to employ the unique properties of 
a frozen mode in a 3D PhC to tailor spontaneous emission. The vanishing group velocity 
enhances emission into the frozen mode, while the perfect impedance match at the PhC/air 
interface leads to efficient extraction of emission from PhC. 

The paper is organized as follows. Section II starts with a brief description of sample 
fabrication process and experimental setup, followed by the experimental results of angle- 
and polarization-resolved reflection and PL spectra. In Section III, we calculate the reduced 
photonic band structure, density of states, and reflectivity to illustrate the effects of frozen 
modes on spontaneous emission process. To show the generality of such effects, we report 
PL enhancement in a different crystal direction in Section IV. Section V is the conclusion. 

II. EXPERIMENT 

A. Sample fabrication 

The fabrication of ZnO inverse opals is detailed in Refs.— Monodisperse polystyrene 
spheres were assembled to face-centered-cubic (FCC) structure on a glass substrate via a 
vertical deposition process. The sample thickness could be varied from 20 to 100 layers 
of spheres. The opal had "domains" with widths of ~ 50 fim and lengths of hundreds of 
micron, separated by cracks. The crystalline arrangement was constant across these cracks 
which formed after the self-assembly process during drying of the opal. ZnO was infiltrated 
into the template by atomic layer deposition (ALD). We ensured the exposure times were 
sufficiently long during ALD growth so that the precursors could fully diffused into the opal 
structure and the ZnO thin film grew conformally and uniformly around each polystyrene 
sphere throughout the sample. The polystyrene spheres were then removed by firing at 
elevated temperature. Figure [T] shows the scanning electron microscope (SEM) images of 
the top surface and cleaved edge of a ZnO inverse opal. The sample surface was parallel to 
(111) crystallographic plane. Even for a 100-layer-thick sample, ZnO was fully infiltrated 
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into the opal template and the filling was nearly 100%. 

A random sample was fabricated to provide reference for optical measurements. 
Polystyrene spheres of different sizes were mixed and deposited on a glass substrate. The 
lack of monodispersity prevented the formation of ordered crystal domains. The randomly 
packed structure was subsequently infiltrated with ZnO and fired at the same temperature 
as for ZnO inverse opal. The identical fabrication conditions ensured similar microstructure 



and material properties. 
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FIG. 1: SEM images of the top surface (a) and cleaved edge (b) of a ZnO inverse opal. 



B. Experimental setup 

We performed the angle- and polarization-resolved refiection measurements to character- 
ize the photonic band structures. Figure [2]^a) shows the experimental setup. The sample 
was mounted on a goniometer stage. The refiection scan was taken in the 6 — 26 geometry, 
namely, when the sample was rotated angle 6, the detection arm was rotated 26. Collimated 
white light from a UV enhanced Xe lamp was focused onto the sample by a lens (LI). The 
beam spot on the sample surface was about 1mm in diameter. The angle of incidence from 
the surface normal was 6. The refiected light was collected by another lens (L2) and focused 
to a fiber bundle which was connected to a spectrometer. A linear polarizer (P) was placed 
in front of the fiber bundle to select s- or p-polarized light with electric field perpendicular 
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FIG. 2: Schematic of experimental setup for angle- and polarization-resolved reflection (a) and PL 
(b) measurements. 

or parallel to the detection plane (made of the detection arm and the normal of sample 
surface). The angular resolution, which was determined mainly by the collection angle of 
lens L2, was about 5°. The reflection spectra were taken with the incident angle 6 varying 
from 5° to 60°. The spectral resolution was about 1 nm. 

In the photoluminescence (PL) experiment, only the detection arm moved and the sample 
did not rotate. As shown in Figure [2](b), the white light was replaced with the He:Cd laser 
light which excited the ZnO. The pump beam was focused onto the sample at a fixed angle 
dp ~ 30°. Spectra of emission into different angle 9 were measured when the detection arm 
was scanned in the horizontal plane. To prevent the reflected pump light from entering the 
detector, the incident beam is shifted vertically so that the incidence plane deviates from 
the detection plane. A long pass filter (F) was placed in front of the fiber bundle to block 
the scattered pump light at wavelength A = 325 nm. 
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C. Reflection spectra 



We measured many ZnO inverse opals with different lattice constants. The sample uni- 
formity was checked carefully with scanning electron microscopy and optical spectroscopy. 
Only the samples that were uniform in crystalline arrangement, thickness, and infiltration 
over an area much larger than the probe beam spot were used in the reflection and PL 
measurements. The data presented below were taken from one sample with the air sphere 
diameter 400 nm and the lattice constant 566 nm. The number of layers of air spheres 
was about 60. Figure [3](a) plots the first Brillouin zone (BZ) of ZnO inverse opal. In the 
reflection and PL measurements, the sample was oriented so that the detector was scanned 
in the TLK plane. Note that in the reflection measurement the UV light from the Xe lamp 
could excite ZnO and generate PL. The PL intensity, however, was much weaker than the 
reflected light intensity. Thus, the PL signals can be ignored in the reflection spectra. 




FIG. 3: (a) The first Brillouin zone (BZ) for a FCC structure. The symmetry points are marked, 
(b) Cross section of the first BZ by the TLK plane. 

Figure ll](a) shows the reflection spectra of p-polarized light for 6 varying from 5° to 55° 
with 5° steps. For 6 = 5°, there is a primary reflection peak centered at wavelength A ~ 800 
nm with the peak reflectivity ~ 0.7. It corresponds to the lowest-order photonic band gap 
in [111] direction, which is confirmed in the numerical simulation to be presented in the next 
section. With increasing angle 6, this peak shifts to shorter wavelength. There are additional 
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reflection peaks at higher frequencies. They are less dispersive and shift to slightly longer 
wavelengths with increasing 9. The reflection spectra of s-polarized light, shown in Fig. 
m^b), exhibited significant difference from those for p-polarized light. There is an additional 
reflection peak at A ^ 510 nm for = 5°. It shifts significantly to longer wavelengths with 
increasing 9. Since the primary reflection peak shifts to the opposite direction, these two 
peaks exhibit an anti-crossing at ~ 45°. 
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FIG. 4: Angle-resolved reflection spectra of a ZnO inverse opal with sphere diameter 400 nm. 
Both the incidence plane and detection plane are parallel to the TLK plane. The incidence angle 
9 ranges from 5° to 50°. The values of 9 are written in the graph. The spectra are offset vertically, 
(a) p-polarized light, (b) s-polarized light. 



D. Photoluminescence data 

When pumped by the He:Cd laser at A = 325 nm, the ZnO inverse opal has PL in both 
the UV and visible frequencies. The UV emission is ascribed to electron transition from the 
conduction band to valence band, and the visible emission is via the defect states within 
the electronic bandgap.— '^'^'^'^'^'^ At room temperature the defect emission, which re- 
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suits from various material defects such as oxygen vacancies, zinc interstitials and oxygen 
interstitials, covers a 300 nm wavelength range. In this paper we concentrate on the defect 
emission, whose broad spectral range allows us to observe the effects of many different-order 
photonic bands on emission. 

Figure shows the PL spectra of the ZnO inverse opal and the random sample. The 
spectral shape of PL from the random sample does not change with observation angle 9. 
The PL of ZnO inverse opal is strongly modified, and the modification is angle- dependent. 
For comparison, the emission spectra taken at identical 9 are scaled so that they overlap at 
A = 880 nm, well below the lowest-order PBG. Suppression of emission at longer wavelength 
is evident in the ZnO inverse opal. At shorter wavelength there is significant enhancement 
of s-polarized emission. This enhancement is not related to stimulated emission, as the 
emission intensity is confirmed to vary linearly with pump intensity. 

We extracted the PL enhancement factor by dividing the PL spectrum of ZnO inverse 
opal by that of the reference sample for same 9. If the normalized PL intensity is less (or 
more) than unity, the spontaneous emission is suppressed (or enhanced). Figure [6] shows 
the normalized PL spectra for both s- and p-polarizations with 9 varying from 0° to 50°. 
The spectra are vertically shifted with a constant offset of 1.5. The reference line of unity 
for each spectrum is plotted as a dashed line. For 9 = 10°, the dip at A ~ 800 nm coincides 
with the primary reflection peak in Fig. |H Its blue-shift with increasing 9 is identical to 
that of the reflection peak. As mentioned earlier, the primary reflection peak corresponds to 
the lowest-order PBG. This partial gap suppresses the emission due to depletion of density 
of states (DOS) within certain angle range. The peaks at higher frequencies in the normal- 
ized PL spectra reveal the emission enhancement by higher-order photonic band structures. 
Most enhancement peaks for both polarizations are weakly dispersive with angle, except 
one for s-polarized PL. This peak red-shifts dramatically with increasing 9. Its amplitude 
reaches a maximal value of 2.3 at ^ = 40°, exceeding all other peaks. It is responsible for 
strong enhancement of s-polarized PL in Fig. [5l We notice that this PL peak has similar 
dispersion to the reflection peak that exists only for s-polarization in Fig. HI To compare 
their frequencies, we overlay the normalized PL spectra and reflection spectra in Fig. [7] (1st 
row) for 9 = 20° — 50°. It is evident that the two peaks do not overlap spectrally, instead 
the PL peak is always at the low frequency shoulder of the reflection peak. 
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FIG. 5: (Color online) Measured PL spectra of a ZnO inverse opal (sphere diameter = 400 nm) 
and a random sample (black dotted line). The emission angle 6 = 20° (a), 30° (b), and 40° (c). 
The red solid curve and blue dashed curve represent the s and p-polarized emission from the ZnO 
inverse opal, respectively. 
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FIG. 6: Normalized PL spectra of the ZnO inverse opal at various emission angles 0. The values 
of 9 are written in the graph. The spectra are vertically shifted with a constant offset of 1.5. The 
reference line of unity for each spectrum is plotted as a horizontal dashed line, (a) p-polarized 
emission, (b) s-polarized emission. 

III. NUMERICAL SIMULATION 
A. Reduced photonic bands 

To interpret the experimental results, we performed numerical simulations of ZnO inverse 
opals using the experimental values. To account for the angular dependence of reflection and 
PL, we calculated photonic bands and density of states (DOS) corresponding to a specified 
angle of incidence/exit in air. These bands and DOS are called reduced bands and reduced 
DOS^-, as opposed to angle-integrated bands and DOS. Since the top surface of ZnO inverse 
opal corresponds to (111) crystallographic plane, we define a reference frame shown in Fig. 
[3l the origin at the F point, z-axis in the TL direction normal to (1 1 1) surface, xy plane 
parallel to the surface with x-axis in the LK direction and y-axis in the LW direction. The 
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wavevector of an incident /exit photon in air can be expressed as 




(sin 6 cos 0, sin 6 sin 0, cos (p) 



(1) 



where uj is the frequency, 6 is the polar angle, (p is the azimuthal angle. It should be 
noted that 9 and denote the direction of light incidence/exit in air, which differs from the 
direction of light propagation inside the PhC Even if the incidence/exit direction in air is 
fixed, the propagation direction inside the PhC may change with frequency uj. 

The reduced photonic band structure is obtained from the full band structure by using two 
conservation laws: (i) parallel momentum conservation modulo a reciprocal lattice vector, 
(ii) energy conservation. They can be written as 



where ky is the k-component parallel to the (111) surface (xy-plane), Gy is the parallel 
component of any reciprocal lattice vector G, and era(k||, k^) represents the nth energy band 
of PhC with wavevector k = ky + k^. 

We used the computer program MULTEM, which is based on the layer KKR (Korringa- 
Kohn-Rostoker) method^, to calculate the reduced band structure. One advantage of this 
method is that it gives k^ for given e„ and ky. Another advantage is that it can calcu- 
late angle-resolved reflection and transmission spectra of a PhC with finite thickness, thus 
allowing direct comparison with the experimental data. The reduced DOS is inversely pro- 
portional to the group velocity, which is calculated from the slope of dispersion curve for 
individual reduced band. Since the frequency range of calculation is below the ZnO elec- 
tronic bandgap, the dispersion of ZnO refractive index is rather weak. We neglected the 
frequency dependence of refractive index in the numerical simulation. To take into account 
the interstitial tetrahedral pores in the ZnO inverse opals^^, we set the value of refractive 
index at 1.95, which is slightly lower than the actual refractive index of ZnO. 

We calculated the reduced photonic bands for kj scanning in the TLK plane, namely 
= and 9 varying from 0° to 80° with 10° steps. The first BZ has mirror symmetry with 
respect to the TLK plane [Fig. [3]^a)]. The s- or p-polarized light, with kj in the TLK plane 
and electric field perpendicular or parallel to the TLK plane, has distinct mirror symmetry 



k|| + G|| = — (sin 9 cos 0, sin 9 sin 0, 0) , 



(2) 



and 



e„(k||,k^) = hu 



(3) 
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with respect to the TLK plane. Hence, the two polarizations are decoupled and the reduced 
photonic bands with kj in the TLK plane are either s- or p-polarized. We plot the s-polarized 
bands in Fig. [7] (2nd row) and p-polarized bands in Fig. [H] (2nd row). The reduced bands 
with positive differ from those with negative k^. This is because the cross section of the 
first BZ by the TLK plane, shown in Fig. [3](b), does not have mirror symmetry with respect 
to X-axis. 

Comparing the reduced s- and p-polarized bands reveals their significant differences. For 
example, at 6 = 30°, the second and third s-polarized bands [labeled as 2s and 3s in Fig. 
[7] exhibit frequency anti-crossing around /c^a/27r = 0.7, while 2p and 3p bands nearly cross 
in Fig. [HI For both polarizations, the second and third reduced bands originate from Bragg 
diffraction of light by (1 1 1) and (1 1 1) planes for k^ > 0, and by (1 1 1) and (2 0) planes 
for kz < 0.^*^ At their crossing point near kza/2Tx ~ 0.7, simultaneous Bragg diffraction by 
(1 1 1) and (1 1 1) planes results in band repulsion.^ The anti-crossing of 2s and 3s bands 
indicates strong band coupling via multiple diffraction. The interaction of 2p and 3p bands, 
however, is much weaker. Such difference can be explained by the dependence of diffraction 
efficiency on polarization. It is well known for X-ray diffraction that when the Bragg angle 
is close to 45° the intensity of diffracted beam is extremely weak for p-polarized wave.^^ The 
suppression of Bragg diffraction has formal analogy to the Brewster effect on reflection of 
p-polarized light by homogeneous medium. The diffraction efficiency for p-polarized light in 
a PhC can be greatly reduced if the incident angle approaches the critical angle.— The 
weak coupling of 2p and 3p bands is attributed to low efficiency of multiple diffraction of 
p-polarized light at the band crossing point because one of the Bragg angles is close to the 
critical angle. 

Comparison of the calculated band structure to the measured reflection spectra con- 
firms that the primary reflection peak corresponds to the lowest-order gap between the first 
and second bands for both polarizations. The additional reflection peak, observed only 
for s-polarization, overlaps with the gap opened by anti-crossing of 2s and 3s bands. As 
9 increases, the gap moves toward lower frequency, and the reflection peak follows. This 
reflection peak is not observed for p-polarization because 2p and 3p bands have little repul- 
sion. In fact, the calculated reflection spectrum for 9 = 40° exhibits a very narrow peak 
corresponding to the small gap between 2p and 3p bands at A;^a/27r ~ 0.7 [Fig. [8]. Such 
narrow peak is smeared out experimentally by averaging over finite collection angle. 
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FIG. 7: (Color online) 1st row: angle-resolved s-polarized reflection spectra (black dashed line) of 
the ZnO inverse opal overlaid with the normalized PL spectra (red solid line) of same polarization 
and angle 9. 2nd row: calculated s-polarized reduced band structure of the ZnO inverse opal for 
a fixed angle of incidence/exit in air. Lattice constant a = 566 nm, the dielectric constant of ZnO 
is e = 3.8. 3rd row: calculated reduced density of s-polarized states of the ZnO inverse opal. 4th 
row: calculated reflectivity of s-polarized light from a ZnO inverse opal whose thickness is 34 layers 
of air spheres. = 0, and 6* = 20° (1st column), 30° (2nd column), 40° (3rd column), 50° (4th 
column). For 6 = 40°, a stationary inflection point is developed for the 2s band at oja/2Tic = 0.856. 

B. Stationary inflection point 

The dispersion of the 2s band is nearly flat in the vicinity of its avoided crossing with the 
3s band. It produces a peak in the reduced DOS shown in Fig. [71 The enhanced PL peak, 
which is observed only for s-polarization, coincides with this DOS peak. It follows the DOS 
peak as it moves to lower frequency at higher 9. Ki 9 = 40°, the dispersion curve for 2s 
band has a stationary inflection point at kzajl-R ~ 0.6 where duj /dkz — and (fu /dkl ~ 0. 



13 



= 20° = 30° = 40° = 50° 




0.7 0.9 1.1 0.7 0.9 1.1 0.7 0.9 1.1 0.7 0.9 1.1 

(0 a/27rc 



FIG. 8: Same as Fig. [7] for p-polarization. 

The existence of stationary inflection point is verified by tracing the evolution of the 2s 
band with 9. For 9 = 50°, the dispersion curve for the 2s band exhibits a local minima at 
kza/27c ~ 0.5 and a local maxima at kza/2'K ~ 0.7. As 9 decreases, the local minima and local 
maxima approach each other, eventually they merge at ^ ~ 40°. With a further decrease of 
9, e.g. at = 30°, the dispersion curve has neither local minima nor local maxima, its slope 
does not change sign throughout the region of interest. Such band evolution confirms that 
the 2s band has a stationary inflection point at = 40° where the merging of a local minima 
and a local maxima gives not only duj/dk^ = but also d?uj/dk1 = 0. The evolution of the 
2p band with 9, shown in the 2nd row of Fig. [8], is completely different from that of the 2s 
band. It reveals that the 2p band does not have a stationary inflection point near 9 = 40°, 
possibly due to its tiny anti-crossing with the 3p band. 

The DOS at the stationary inflection point diverges in an inflnite large PhC In a real 
sample such divergence is avoided because of finite sample size. Nevertheless, the DOS peak 
has the maximal amplitude at ^ ~ 40° where the stationary inflection point is developed. 
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The large DOS enhances the spontaneous emission process. Although the 2p band does 
not have a stationary inflection point, its dispersion is relatively flat in the neighborhood 
of kza/2TT = 0.6 [Fig. [S]. It produces a peak in the reduced DOS, which should enhance 
emission. Experimentally, the p-polarized PL is not enhanced. 

The question arises why the DOS peak leads to enhanced PL for s-polarization but not for 
p-polarization. The answer lies in the emission extraction efficiency. One unique property 
of the frozen mode at the stationary inflection point is its efficient coupling to the free 
photon mode outside the PhC It leads to vanishing reflectivity at the sample/air interface, 
which is conflrmed by our calculation and measurement of reflectivity from ZnO inverse 
opal. The calculated reflection spectrum for 6 = 40° [Fig. [7] reveals that the reflectivity at 
the stationary inflection point uja/2'Kc = 0.856 is almost zero. Experimentally, the measured 
reflectivity is low but not zero due to averaging over flnite collection angle. To verify that the 
zero reflectivity is not a result of flnite sample thickness, we calculated the reflection spectra 
of three ZnO inverse opals with different thicknesses (34, 36 and 38 layers of air spheres). 
Figure M^a) shows that for 6 = 40° the reflectivity reaches zero at multiple frequencies. 
As the sample thickness varies, all the zero points of reflectivity shift in frequency except 
the one at uja/2'Kc = 0.856. Their dependence on sample thickness suggests those vanishing 
reflectivities result from interference of light multiply reflected by the two surfaces of the ZnO 
inverse opal. The fact that the zero reflectivity at uja/2'Kc = 0.856 is independent of sample 
thickness conflrms that it is caused not by the Fabry-Perot resonance in a flnite PhC slab 
but by the intrinsic property of photonic band, more specifically, the dispersion of 2s band. 
This result demonstrates a perfect impedance match for the frozen mode at the interface of 
ZnO inverse opal and air. The light emitted to the frozen mode experiences little reflection 
at the sample surface when leaving the sample. Hence, s-polarized PL at the stationary 
inflection point of 2s band can be efficiently extracted from the sample. The reflection of 
]9-polarized light does not go to zero in the absence of frozen mode, thus p-polarized PL 
cannot escape easily from the sample. 

Therefore, the observed large enhancement peak of s-polarized PL around 6 = 40° is 
attributed to two factors: (i) enhanced emission into the frozen mode due to large DOS, 
(ii) efficient extraction of emitted light out of the PhC. As 6 approaches 40°, the DOS 
increases and surface reflection decreases. Thus, PL is enhanced in the vicinity of a stationary 
inflection point. The maximal PL enhancement occurs at the stationary inflection point 
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FIG. 9: (Color online) (a) Calculated reflectivity of s-polarized light as a function of normalized 
frequency ujaj^ixc. The incidence angle Q = 40°, cj) = 0. The thickness of ZnO inverse opal is equal 
to 34 (red solid line), 36 (blue dashed line) and 38 (black dotted line) layers of air spheres, (b) 
Calculated reflectivity of s-polarized light at the stationary inflection point uja/2TTC = 0.856 as a 
function of the azimuthal angle (p. 9 is fixed at 40°. 

where the DOS reaches the maximum and the surface reflection the minimum. 

The enhanced emission is directional, namely, it exits the ZnO inverse opal to air at the 
polar angle 6 = 40°. Next we investigate the emission directionality in terms of azimuthal 
angle 0. Since the cross section of the first BZ by the TLK plane does not have mirror 
symmetry with respect to the z-axis (FL direction), the reduced band structure for kj 
scanning along the LK path is different from that along the LU path. The reflection spectra, 
however, are the same for the two scanning directions due to reciprocity of reflection.— This 
is confirmed by our calculation and measurement of reflection spectra. Figure [9](b) plots 
the calculated reflectivity for s-polarized light at the frequency of stationary inflection point 
versus the azimuthal angle 0. The reflectivity exhibits six-fold symmetry when varies from 
0° to 360°. In principle, the PL spectra should exhibit difference between the scanning path 
along LK and that along LU. Experimentally, the PL spectra for the two scanning paths 
are similar. This is attributed to the structure disorder, in particular, the twin structure 
in the ZnO inverse opal.^^ As a result, the enhanced PL can be observed at six azimuthal 
angles. 
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IV. PL ENHANCEMENT IN A DIFFERENT CRYSTAL DIRECTION 
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FIG. 10: Measured angle-resolved spectra of normalized PL (a) and reflection (b) of a ZnO inverse 
opal with sphere diameter = 400 nm. The scanning plane is parallel to the TLW plane, (j) = 90°. 
9 varies from 10° to 60°. The values of 6 are written on the graph. In (a), the normalized PL 
spectra are shifted vertically with a constant offset of 2. The horizontal dashed lines mark unity 
for individual spectra. The reflection spectra in (b) are also shifted vertically. 

Although the data presented above were taken from the ZnO inverse opal with sphere 
diameter 400 nm, we repeated the experiments with several samples of different sphere size 
and obtained similar results. The enhancement of PL by the frozen mode is a common 
phenomenon because many high-order bands of ZnO inverse opals have stationary inflection 
points in their dispersion curves. In addition to scanning in the TLK plane, we also scanned 
in the TLW plane and observed PL enhancement at a different angle 6. 

Fig. [TOT a) shows the angle-resolved PL spectra when scanning in the TLW plane. The 
normalized PL spectra are shifted vertically with a constant offset of 2. We observed an 
enhancement peak which shifts to lower frequency with increasing 6. The maximum en- 
hancement occurred around 9 = 50° and reached the value of 3.2. Since the cross section of 
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FIG. 11: (Color online) 1st row: angle-resolved reflection spectra (black dashed line) of the ZnO 
inverse opal overlaid with the normalized PL spectra (red solid line) of same angle 9. 2nd row: 
calculated reduced band structure of the ZnO inverse opal for a fixed angle of incidence/exit in air. 
Lattice constant a = 566 nm, the dielectric constant of ZnO e = 3.8. 3rd row: calculated reduced 
density of states of the ZnO inverse opal. 4th row: calculated reflectivity of a ZnO inverse opal 
with thickness = 34 layers of air spheres, (j) = 90°, and 9 = 30° (1st column), 40° (2nd column), 
50° (3rd column), 60° (4th column). For 9 = 50°, a stationary inflection point is developed at 
uja/27:c = 0.87. 

the first BZ by the TLW plane is symmetric with respect to y-axis (parallel to FL direction), 
the photonic mode with (k|[, k^) is identical to that with (k||, — A;^). The reduced band struc- 
ture with positive kz is the same as that with negative kz- However, the first BZ does not 
have mirror symmetry with respect to the TLW plane. The photonic modes with kj in the 
TLW plane contain both electric field components parallel and perpendicular to the TLW 
plane. Hence, the reduced photonic bands can no longer be divided into s- and p-polarized 
bands. Consequently, the PL spectra as well as the reflection spectra are insensitive to 



18 



polarization. The angle-resolved reflection spectra in Fig. [TO](b) exhibit the frequency shift 
of the first and second refiection peaks with 9. Comparing the refiection spectra to the PL 
spectra in Fig. [TT] reveals that the enhanced PL peak overlaps with the refiection dip at 
the low-frequency side of the second refiection peak. Our calculation of the reduced band 
structure, the reduced DOS and refiection spectra (Fig. [TTi) illustrates that the maximal PL 
enhancement peak observed around 9 = 50° results from the stationary infiection point of a 
high-order band, which produces a peak in DOS and a dip in refiectivity at u;a/27rc ~ 0.87. 
Note that the refiectivity at 9 = 50° is not zero, because the exact angle for the stationary 
infiection point is 47°. The comparison of numerical results to experimental data confirm 
that the PL enhancement originates from the frozen modes, similar to that observed when 
scanning in the TLK . 



V. CONCLUSION 



In conclusion, we present a detailed study on the angle- and polarization-resolved pho- 
toluminescence and refiection spectra of ZnO inverse opals. The broad ZnO defect emission 
exhibits multiple enhancement peaks for both polarizations. Our numerical simulation re- 
veals that the largest enhancement peak results from the frozen mode at the stationary 
infiection point of dispersion curve for a high-order photonic band. The frozen mode has 
well defined propagation direction and may be linearly polarized. At its frequency, the 
reduced DOS diverges, greatly enhancing the spontaneous emission into the frozen mode. 
Perfect coupling of the frozen mode to the free photon mode outside the sample leads to 
efficient extraction of emission from the PhC The enhanced emission not only has good 
directionality but also can be linearly polarized. We note that the above mechanism does 
not work for all the enhanced PL peaks. Many peaks, which overlap spectrally with the 
high-order photonic bands, remain to be explained. Nevertheless, our results demonstrate 
that a high-order band of a 3D PhC can strongly modify the spontaneous emission, which 
offers potential application for a highly efficient light source. 
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